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The evaporation of sessile drops in quiescent air is usually governed by vapour diffusion. 
For contact angles below 90°, the evaporative flux from the droplet tends to diverge in 
the vicinity of the contact line. Therefore, the description of the flow inside an evaporat- 
ing drop has remained a challenge. Here, we focus on the asymptotic behaviour near the 
pinned contact line, by analytically solving the Stokes equations in a wedge geometry of 
arbitrary contact angle. The flow field is described by similarity solutions, with exponents 
that match the singular boundary condition due to evaporation. We demonstrate that 
there are three contributions to the flow in a wedge: the evaporative flux, the downward 
motion of the liquid-air interface and the eigenmode solution which fulfils the homo- 
geneous boundary conditions. Below a critical contact angle of 133.4°, the evaporative 
flux solution will dominate, while above this angle the eigenmode solution dominates. 
We demonstrate that for small contact angles, the velocity field is very accurately de- 
scribed by the lubrication approximation. For larger contact angles, the flow separates 
into regions where the flow is reversing towards the drop centre. 



1. Introduction 

Evaporation of colloidal dispersion droplets is a widely-used mechanism to deposit 



particles onto a substrate and generate colloidal crystals (Bigioni et al. 2006 Dufresne 



et aL|2003 Velikov|2002 1 . The ring-shaped stains that remain after evaporation can also 



be disadvantageous, for example in the coating and inkjet-printing industry (Deegan et al. 



1997 2000 Eral et al. 2011 Yunker et al. 2011 ). To understand and control the stains that 



form when a droplet evaporates, one needs to know the velocity field inside the drying 
drop peegan et Z1|1997| |Marm et q/."|2011| |Ristenpart et a/.||2007| |Brutin et al\2Qll\ . 
The capillary flow inside an evaporating drop is driven by the evaporative mass loss from 
its surface. There are three mechanisms which can be rate limiting for the evaporation 
of a drop (IHaut & Colinet||2005| ICazabat & Guena||2010| IMurisic & Kondic||2011|): the 



transfer of molecules across liquid-air interface, the heat transfer to the interface, or the 
diffusive transport of the vapour in air. One of the first two mechanisms can be dominant 



when thin films of evaporating liquid are considered (Cazabat & Guena 2010 Haut & 



1988 



1995 



Colinet||2005 ), when the surrounding phase is not gas but pure vapour (Burelbach et al. 



Colinet k. Rednikov|2011 ), or for droplets on heated substrates (Anderson & Davis 



Cazabat & Guena 2010 ). For macroscopic evaporating drops in air, diffusion- limited 



evaporation is often assumed, based on estimates of the time scales for transport ( Hu & 



Larson 2002 Popov 2005 Eggers & Pismen 2010 Cazabat fc Guena 2010). For water 



drops there has been some debate about the time scale for transport across the interface, 



and hence about the applicability of diffusion- limited models (Cazabat & Guena 2010 
Murisic fc Kondic|2011[ ). Experimentally however, the diffusion-based evaporation model 
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is found to describe the evolution of the droplet mass and contact angle of sessile water 
drops with pinned contact lines very well, for the entire range of possible contact angles 



(Deegan et al. 


1997 


Hu & Larson 


2002 


Guena et al. 


2007 


Cazabat & Guena 


2010 


Gelderblom et al.\\2011 


Sobac & Brutin||2011 





Here, we will study the flow near the pinned contact line of a macroscopic evaporating 
drop in air on an unheated substrate, and therefore we consider diffusion-limited evapo- 
ration. Until now, the nature of the flow in the vicinity of the contact line has remained 
unclear. In the diffusion-limited case, the singular corner geometry of the droplet close to 
the contact line gives rise to a diverging evaporative flux and hence to a diverging velocity 
field peegan et aZ.||1997[ [2000| |Hu fc Larson||2005| |Popov||2005[ ). This singularity makes 
analytical and numerical solutions to the velocity field inside a drop difficult to obtain 



(Colinet & Rednikov 2011 


Hu & Larson|2005 


Fischer|2002 


Masoud & Felskc 2009 


Petsi 


& Burganos||2008[ Poulard et al. 


2005 


). In several studies the flow inside the drop was 



solved analytically, however, at the expense of smoothing the evaporative flux singular- 



ity. Masoud & Felske ( 2009 ) considered an exponential cut-off for the flux, while Petsi 
& Burganos (2008[) focussed on uniform evaporation profiles. For small contact angles, 



evaporation-driven flow inside a droplet is often described in the lubrication approxi- 



mation (Berteloot et al. 2008 Eggers & Pismen 2010), which compares very well with 



experimental data (Marin et al. 2011). It was argued by Hu & Larson (2005) however, 
that the standard lubrication approximation does not hold near the contact-line region 
due to the diverging evaporative flux. 

On top of that, Marangoni stresses could alter the velocity field in the vicinity of the 
contact line: the non-uniform evaporative flux leads to temperature gradients over the 
drop surface, which give rise to differences in surface tension, and drive a Marangoni 



flow inside the drop, as has been confirmed experimentally by Hu & Larson (2006) for 
octane droplets. Dimensional analysis shows that this Marangoni effect is so strong that 
it can overcome the diverging evaporation-driven outward flow, at least at some distance 
from the contact line ( Bodiguel fc Leng|2010 Hu fc Larson|2006 Ristenpart et a/.|2007 ). 
However, for water droplets the Marangoni effect is found to be weak (Hu & Larson 



2006). In PIV measurements of the velocity field in an evaporating drop by Marin et al. 
(2011) the experimental velocities were of the order of 10 /xm/s, whereas the Marangoni 
velocities would be of the order of 10 mm/s. 

Here, we derive analytical solutions of evaporation-driven Stokes flow in a wedge ge- 
ometry to address the nature of flow near the pinned contact line of an evaporating drop. 
While solutions to the full flow pattern in the drop can only be obtained numerically ( Hu 



& Larson 2005 Fischer 



2008 Masoud & Felske 



2002) or for a regularized evaporative flux (Petsi & Burganos 



2009). the behaviour in the vicinity of the contact line is char- 



acterized by similarity solutions. This is a classical approach for flows near contact lines 



that goes back to Huh & Scriven ( 1971 ) and was recently applied to Marangoni flow in 



evaporating drops by Ristenpart et al. (2007). We examine the velocity field inside the 



drop while retaining the singular evaporative flux as a boundary condition. We demon- 
strate that for small enough contact angle, the lubrication approximation can be applied 



all the way down to the contact line, which invalidates the argument of Hu & Larson 



(2006). For larger contact angles (above 127°) interesting flow structures appear, with a 
reversal in the flow direction. We show that there are three contributions to the total flow 
in the wedge: one that comes from the evaporative flux boundary condition, one from 
the downward movement of the liquid-air interface, and one eigenmode solution, which 
satisfies the homogeneous boundary conditions. Which of these conditions is dominant, 
depends on the contact angle 9 of the drop. For 8 < 6 C = 133.4°, the critical angle the 
evaporative flux solution dominates, whereas for 9 > 9 C the eigenmode solution domi- 
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Figure 1 . (a) An evaporating drop on a substrate with a contact angle 9 and base radius R. The 
arrows indicate the evaporative flux of water vapour J from the drop surface to the surroundings. 
The dashed square marks the area close to the contact line, where the drop geometry can be 
approximated by a wedge, (b) Overview of the wedge geometry. The evaporative flux drives the 
Stokes flow inside the liquid. The contact line is located at the origin of the polar coordinate 
system (p, (f>). 



nates. The solution at the critical point is treated separately. Finally, we comment on 
the typical pressure in the vicinity of the contact line and on the regularization of the 
evaporative singularity. 



2. Corner solutions 

2.1. Problem formulation 

The geometry of the droplet close to the contact line can be approximated by a two- 
dimensional wedge with contact angle 9; see figure [l] The contact line is located at the 
origin of the polar coordinate system (p, 4>) . 

We define the velocities in terms of streamfunction ^(p, </>) as 



(2.1) 



19* , ,. 9* 

The flow is governed by the Stokes equations, or equivalently, in terms of the stream- 
function, by the biharmonic equation 



V 4 * = 0. 

As boundary conditions we have no slip and impermeability of the substrate 
which the droplet is deposited 



u p (p,0) 



i a* 

p d(j) 



= 0, and u,p(p, 0) 



!>=0 



and no stress at the liquid-air interface 



a* 

dp 



= 0, 



/>=o 



(2.2) 
= 0) on 

(2.3) 



(>-. 



8 



i a* 

p dp 



i a 2 * 



o, 



(2.4) 



dp \p op j p* 

with 7/ the dynamic viscosity. The problem is closed by the kinematic boundary condition 
at the liquid-air interface, where mass transfer due to evaporation occurs. This kinematic 
boundary condition consists of two contributions. One contribution comes from the evap- 
orative mass flux from the interface, which drives a flow inside the droplet. Due to the 
evaporative mass loss, the droplet volume decreases with time. While the contact area of 
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the drop remains constant because the contact line is pinned, the contact angle decreases 
with time. This gives rise to a second contribution to the kinematic boundary condition: 
the downward motion of the liquid-air interface acts like a closing hinge. Hence, the 
kinematic boundary condition reads 



dp 



= -J(P) + 
Pi 



dO 
dt 1 



(2.5) 



with J the evaporative flux and pi the liquid density. Both contributions in (2.5) are 
known in detail from earlier studies ( Deegan et aZ.|1997 Popov|20 05). The key ingredient 
is that the mass loss from the droplet is limited by diffusive transport of the water vapour 
in the air outside the drop. By solving the vapour concentration field outside the droplet, 
one can find an expr ession for the evapo rative flux J, which, close to the contact line, 



scales as J - 
and 



_1 ( Deegan et aL|l997 ), where p = p/R, with R the drop base radius, 



MO) 



(2.6) 



2tt-28 

Hence, for A < 1, which means < 90°, the evaporative flux diverges as the contact 
line is approached. From the solution for the vapour concentration field in the corner 
geometry, the evaporative flux is found to be ( Deeg an et aZ.||1997| ) 

J(P) 



Pi 



A{8)Up m ~ 



(2.7) 



where U = DAc/Rpi is the velocity scale, which is of order pm/s for water drops in 
ambient conditions ( |Marm et al. 2011[ ), with D the diffusion constant for vapour in 
air, and Ac = c s — Coo the vapour concentration difference (in kg/m 3 ) between the drop 
surface and the surroundings. Prefactor A{9) can be found from the asymptotic behaviour 
of the full spherical-cap solution ( Popov||2005|) and is of order unity. The rate of contact 



angle decrease, d9/dt, can be determined from the total rate of mass loss from the drop. 
A closed-form analytical solution for the rate of contact angle decrease was derived by 



Popov (2005) 



(W 
~dt 



-B{6) 



U 
i?' 



with 



B(6) = (1 + cos( 



smt 



1 + cosh 29t 



tanh(7r — #)rdr 



(2.8) 



(2.9) 



1 + cos 6 ' J Q sinh 2itt 
2.2. Solution 

For a given power n, a general solution to biharmonic equation §H§ reads ( |Michell|189"9| 
\l/(yO, 4>) = p n [ci cosn</> + C2 sinn(j) + C3 cos(n — 2)(f> + C4 sin(n — 2)</>] , (2.10) 



where n ^ 0, 1,2. Since (2.2) is a linear equation, the two contributions to the inhomoge- 



neous boundary condition (2.5) can be considered separately, whereas the full solution is 
obtained by superposition. On top of this, one may in principle add the classical "eigen- 



mode" solution of the homogeneous problem by Dean & Montagnon (19491; Moffatt 



(1964), for which all boundary conditions are zero. Hence, we need to consider three 



types of solutions: the velocity field due to the evaporative flux condition 



u<p(p,0) = — J(p), 
Pi 



(2.11) 
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, . d6 
u${p,V) = -rrp, 

and the flow that satisfies the homogeneous condition 



(2.12) 



u 4> (p,9) = Q. (2.13) 
From now on, we will refer to (2.11) as the flux condition, to (2.12| as the hinge 



condition, while the solution satisfying (2.13) is the corner eigenmode. Note that each of 



these boundary co nditio ns will give rise to a different p ower of p in the final solution. 
The flux condition (2.11 ) scales as p x ~ 1 , with A given by (2.6), while the hinge condition 



(2.12) scales as p. The homogeneous eigenmode solutions turn out to scale with yet 
another exponent, denoted A^- This exponent also depends on 9, and follows from an 
eigenvalue equation M{X E ,9) = (Moffatt 1964), where 

(2.14) 



M(A, 9) = sin 2(A - 1)9 - (A - 1) sin 29. 



This equation has a countable infinite number of zeros, each of which corresponds to a 
different eigenmode (with the exception of the trivial solutions A = 0, 1, 2). The solutions 



of (|2.14[) are discussed in great detail by |Dean fc Montagnon| ([1949J; |Moffatt| ( |1964[ ); 
Moffatt fe Duffy| ( |1980[ ). For angles 9 < 79.6° (|2.14[ ) has complex roots, which causes 
viscous eddies to appear in the flow QMoffatt||1964[ ). We are interested in the lowest root 
that has TZ(Xe) > 1, to ensure regularity of the eigenmode velocity. 



2.2.1. Flux condition 



One can cast the solutions generated by the flux condition (2.11 ) in the form 



RUA{9) . 



P x f{^0), 



M{\,9) 

where A depends on 9 as given in (|2.6|), and the ^-dependent part reads 



(2.15) 



/(& 0) = 2 {( A ~ 2 ) t sin Xe ~ sin ( A _ 2 ) 6 '] t cos X< t> ~ cos ( A _ 2 )$ 



+ [A cos X9 - (A - 2) cos(A - 2)0] 



A - 2 

sin(A — 2)<p — sin A</> 

A 



(2.K 



The factor RU provides the dimensional strength of the streamfunction, while A(6) 
captures the dependence of the evaporative flux on the contact angle of the drop. In- 
terestingly, the denominator contains the factor M(A, 9) that was previously defined in 
(2.14). We thus need to separately consider the cases where M{\,9) = 0, for which the 
solution (2.15) is not defined. 



The function M(\, 9) has two obvious roots that are encountered in the flux problem, 
namely A = 1 and A = 2. These correspond to cases for which the form (2.10) is degen- 
erate and additional solutions to the biharmonic equation appear. For A = 1 (6 = 90°) 
the flux solution becomes 



j * 2RUA(w/2) _ . 
*(p, 4>) = v / ' p&smt 



(2.17) 



while for A = 2 (9 = 135°) we find 
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M/A(3tt/4) 



~p z (1 - cos 2</>) . 



(2.18) 



90° and 



Note that these degenerate solutions are a regular limit of (2.15) for 9 
135°, respectively. 

In addition to these trivial roots, M(A, 6) with A = tt/2(tt — 9) exhibits one root 
that leads to a truly nontrivial solution. This root appears at a critical angle 9 C = 133.4°, 
with corresponding exponent A c = 1.93. The critical point arises when the four boundary 
conditions are not linearly independent, which implies that the inhomogeneous system 
cannot be solved. Indeed, this is exactly the condition required for a nontrivial (eigen- 
mode) solution of the homogeneous problem, namely M(Xe,9) — 0. As a consequence, 
the eigenmode solution has the same exponent A# = A c at the critical angle 9 C . The 



resulting critical solution is not of the form (2.10) and will be treated separately in 5 2.3 



2.2.2. Hinge condition 



The hinge condition (2.12) gives rise to a degeneracy of solution (2.10), since in this 



case n = 2. The solution to (2.2) with the hinge condition is given by 



T . ,. RUB{9) _ 2 



N(0) 



with N(9) = 2(29- tan 20), and 



9( 



= sin 2(j> - tan 29 cos 20 - 2(f) + tan 29. 



(2.19) 



(2.20) 



Hence, (2.19) is defined for all angles, except at 9 — and 9 = 9h = 128.7, where 



N(9) = 0. For 9 satisfying N(9) = 0, A = 2 becomes a double root of eigenvalue equation 
(2.14). Hence, once more a critical point appears when the eigenmode has the same 
exponent as the solution of the inhomogeneous hinge problem, i.e. A^ = 2. We anticipate 
that this critical point will be of less importance than the critical point for the flux 
condition, as the latter has a slightly smaller exponent A < 2. 



2.2.3. Eigenmode 

The eigenmode represents the nontrivial solution of the homogeneous problem, ap- 
pearing when M(\ E ,0)= 0. It reads ( |Moffatt||1964[ ) 



9(p, <t>) = RUC{6)p x *h{<l>,e), 



with 



h{<t>,e)= S in[{\ E -2){e-cj>)]- 



sin(A i5 - 2)6* 



sm[\ B {e - cj>)]. 



(2.21) 



(2.22) 



sin Ab£ 

The prefactor C(9) cannot be determined from the "inner" Stokes flow problem, where 
we consider only the wedge in vicinity of the contact line. This is fundamentally different 
from the prefactors A(0), B(9), which are known from the external boundary conditions 
on the wedge. In general, the eigenmode solution will be excited by the far-field flow 
inside the drop, and is therefore determined on the outer scale R. This is beyond the 
present, local analysis. As the only velocity scale in the problem is the one induced by 
the evaporation, the streamfunction will naturally scale as RU, and C(9) will be of order 



unity (with the exception of the critical point). Note that (2.21) is the classical solution 



by Moffatt ( 1964 ) , which leads to the famous viscous eddies when A e has a nonzero 



imaginary part. 
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2.3. The critical point 

The vanishing denominator of ( |2.15| ) at 9 C = 133.4° and ( |2.19[ ) at 6 h = 128.7° signals 
the breakdown of the local similarity solution. Such a breakdown of similarity solutions 
in corner flows has been analysed in great detail by Moffatt & Duffy (1980). This work 
considered a pressure-driven flow along a duct which cross-section has a sharp corner, as 
well as a variation to the hinge problem considered above. All cases displayed the same 
scenario: the inhomogeneous boundary conditions cannot be fulfilled at a critical angle 9 C , 
due to an overlap with the homogeneous eigenmode. The key result by Moffatt fc Duffy| 
(1980) is that the critical solution develops logarithmic corrections, of the type ^f c ~ 
p Ac In p. This can be derived by considering the regular solution, including the eigenmode 
contribution, in the limit 9 — > 9 C . We therefore introduce an expansion parameter, e = 
9 — 6 C , and derive the critical flux solution in the limit e — > 0. Below we discuss in 
detail the critical point for the flux condition, as it will turn out the most relevant 
for the evaporation problem. The critical point for the hinge condition can be treated 
analogously; the result is given below. 

The identity of exponents at the critical point, A = A^ = A c , forces us to consider a 
superposition of the flux solution and the eigenmode, f c — ^ f + ^e- As expected, the 
flux solution diverges near the critical point as 1/e, and gives an expansion of (2.15) 



^ f = ^^[A + eA 1 + O(e 2 )] [/ o (0) + e/i(0) + O(e 2 )]. (2.23) 

Here we note that A = A c + eAi + O (e 2 ), and all coefficients can in principle be de- 
rived from the expressions given in this section. Here we summarize the leading order 
contributions 



2(7T-0 C )2> 



(2.24) 



where f(4>,9) was previously defined in ( |2.16 1, and Mx(X c ,0 c ) — dM/d9\g c . As noted by 
Moffatt & Duffy (1980), a regular solution for e — > is only achieved if the eigenmode 



displays an identical 1/e scaling, to compensate for the divergence of tyf. To leading 
order, we thus require C(ff) ~ —A /e, such that the expansion of the eigenmode can be 
written as 



RUp> 



[Ao +eC\+0 (e 2 )] [ho(<f>) + ehitf) + O (e 2 )] , 



(2.25) 



where now A# = A c + eA^i + O (e 2 ). Once again, we provide the leading order contri- 
butions 



h Q {4>) = h(<p,9 c ), X E ,i = 



~a¥ 



2(A C - 1) [cos 2(A C - 1)6 c - cos 29 c 
29 c cos 2(A C - 1)6 c - sin 26» c 



(2.26) 



Indeed, one can verify that ho(<p) = fo(4>), which ensures a perfect cancellation of the 
1/e contributions of and ^e- Finally, the critical solution *1> C = ^ f + is obtained 
from ((27231) and (p725| as 
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^ c = RUp x " 



x <^ A 



/„(0) + (A! - Ci)/ o (0) + A (A(0) - fti (</»)) + O (e) 



#[M 0/ 5 A = {(Ai - X Etl ) In (p/«) /„(0 + [AW - + (e)} . 



(2.27) 



Indeed, this confirms the scenario that the leading order asymptotics of the critical 
solution is of the form ^f c ~ p x " \n(p / k) . The logarithm appears due the expansion of 
~eAi _ pcA E .i^ an( j tnug or igi na tes from the "closeness" of A and A^ near the critical 
point. The length scale k that appears inside the logarithm cannot be determined from 
the current local analysis. Namely, k follows from the combination (Ai — C\) appearing 



in (2.271. The coefficient C± requires more knowledge of the eigenmode amplitude C{9), 



and thus of the large-scale flow in the spherical-cap shaped drop. Finally, one can verify 
that this critical solution indeed satisfies the inhomogeneous boundary condition, due to 
the properties A fi — A(9 C )/X C and /q = hx = at the free surface = 9 C . 

The hinge solution at the critical point 0^ = 128.7° can be obtained following a similar 
procedure: the singular contributions from the hinge and eigenmode solutions at this 
point will give rise to logarithmic corrections. Here we only state the result 



RUB{9 h ) _ 2 

o/l3 P 



In (p/K h )g ((p) 



[ffi^-AiMl+Ol 



(2.28) 



with go((f>) = g{4>,9h) = sin2</> — 2cf> + 20^(1 — cos2</>). Note that the leading order part, 
which contains the logarithmic corrections, is identical to the hinge problem considered 



by Moffatt & Duffy (1980). Differences arise in the inhomogeneous contribution, gi{4>), 
due to the different boundary conditions. 



3. Results 

3.1. Dominant contribution 

The complete flow field is obtained by a superposition of the flux, hinge and eigenmode 
solutions identified above. Which of these will be relevant near the contact line depends 
on the scaling with p: the lowest exponent provides the leading order asymptotic con- 
tribution. Figure [2] shows the exponents A (flux solution, solid line) and Xe (eigenmode, 
dotted line), as a function of the contact angle 9. For 9 > 9 C the flux solution has the 
higher exponent, and hence the eigenmode will provide the leading order contribution. 
For 9 < 9 C) the flux solution dominates. The hinge solution ~ p 2 (dashed line) is asymp- 
totically subdominant for all contact angles. Note that in the vicinity of 9 C — 133.4°, 
however, the values of all exponents are very close to 2. This means that in this range of 
contact angles the asymptotic solution can be reached only at very small p: on practical 
scales all solutions will contribute. 

Most experiments on evaporating drops with pinned contact lines are performed at 
angles below 90°. This means that the flux solution, for which the prefactor A(9) is 
known, is by far the most relevant case. Nevertheless, an important conclusion is that a 
local analysis of the problem cannot provide the amplitude of the leading order flow for 
9 > 9 C : the prefactor C(9) is determined from matching to the outer flow. 

3.2. Streamlines 



The solutions (2.15), (2.19) and (2.21) in principle contain all information on the liquid 
flow in the vicinity of the contact line. We plot the streamlines associated to these so- 
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100 
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Figure 2. A plot of the exponents of p that arise in the flux solution (solid line), hinge solution 
(dashed line) and eigenmode solution (dotted line) versus 9. For 9 < f33.4°, the flux solution is 
dominant, whereas for larger 9 the eigenmode solution dominates. The hinge solution is always 
subdominant. The critical points 9 C = 133.4 and 9h = 128.7 arise when respectively the flux or 
hinge solution is equal to the eigenmode solution. 




Figure 3. Streamline plot of the flux solution ( |2.15[ ) for different contac t an gles. For 9 — 128°, 
130°, 150° and 160° reversal of the flow direction is observed. Based on ( |3.2[ ), a new separatrix 
(dashed line) appears for 9 m 127° (which disappears at 9 = 9 C — 133.4°), 9 R3 148° and 
9 « 157°. At 9 — 9 C we only show the (^-depen dent part of the solution, since the prefactor is 
diverging and has to be treated separately (see j|2.3[). 



lutions in figure [3j [4] and [5] for different values of the contact angle. Figure [3] shows the 
flow due to the flux condition, with a mass transfer out-of the liquid, figure [4] represents 
the hinge condition due to the moving liquid-air interface, while figure [5] shows the eigen- 
mode solutions. For the inhomogeneous solutions, the streamlines arrive at the liquid-air 
interface with a well-defined angle. As we deal with similarity solutions, the inclination of 
the streamlines with respect to the free surface is independent of p and depends only on 
9. The only exception is the solution at the critical point 6 C (2.27), where an additional 



length scale k from the outer problem comes in and the self-similarity is lost. 

When analysing the flux solution in figure[3]in more detail, one can see that for < 90° 
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Figure 4. Streamline plot of the hinge solution (2.191 for different contact angles. A separatrix 
(dashed line) appears for 9 = 90°, and disappears again for 9 = 9h ~ 128.7°. At 9 = 9h we only 
show the ^-dependent part of the solution, since the prefactor is diverging and has to be treated 
separately (see \ 2.3 1 



6=30" 




8=60" 




9=90" 




Figure 5. Streamline plot of the eigenmode solution (2.211 for different contact angles with 
prefactor C{6) — 1. For 6 = 30° and 60° viscous eddies are present. The dashed line separates 
the two eddies visible for 9 — 30°. 



the strength of the flow increases upon approaching the contact line. From the scaling 
in (2.151 one in fact sees that the velocity diverges as p — » for this range of contact 
angles, since A — 1 < 0. Above 8 > 90°, however, the velocity decays and vanishes near 
the contact line. Interestingly, the flow displays some reversal structure for these large 
contact angles. As 9 increases from 126° to 128° the flow becomes separated into two 
regions, such that near the bottom wall the liquid flow direction is actually away from 
the contact line, towards the center of the drop. The dashed line shows the separating 
flow line that ends with a stagnation point at the contact line. The separatrix then moves 
upwards with increasing 9, as seen in the plot for 9 = 130°. Between 9 = 130° and 140° 
the separatrix disappears. This occurs when the separatrix reaches the free surface, which 
coincides with the critical point, i.e. 9 = 9 C = 133.4°. At this critical angle, the surface 
flow changes direction. Indeed, a separatrix located at the free surface is incompatible 
with the boundary condition of an outward flux, which illustrates the breakdown of a 
local similarity solution and the appearance of the critical solution (2.27). When the 
contact angle is increased to 150° a new separatrix appears. For 9 — 160° we see this 
separatrix has moved upwards, and a second separatrix has appeared. Separatrices that 



appear for 9 > 133.4° do not disappear again, since the flux solution (2.15) has only 



one critical point. One can demonstrate from the exact solutions that a new separatrix 



appears at 4> — when du p /d(f>\ 



<P=o 



0. This leads to the condition 



p{9) = sin(9(A - 2) - sinflA = 0, 



(3.1) 



with A = 7r/2(7r — 9), which is illustrated in figure [6] From criterion (3.1) the contact 
angles at which a new separatrix appears are found to be (in radians) 



1 + 6fc + k 2 + 1 - k , for k = 1, 3, 5 



(3.2) 
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Figure 6. The contact angles at which a new separatrix appears for the flux condition are 
located at the zeros of p(9) = sin \8 — sin(A — 2)9, given by ( 3.2 1. 



The hinge solutions have a simpler structure (figure [2]). For all values of 9 we find that 
u r*> p, and thus the speed increases linearly with the distance from the contact line. For 
contact angles below 90° the flow is oriented away from the contact line. For 9 > 90° only 
one single separatrix appears, which has disappeared again for 6 > 128.7°. At the critical 
point Oh = 128.7° the local similarity solution again breaks down, and we encounter a 
critical solution. 

The eigenmode solutions are shown in figure [5j where we took C{9) = 1. In the 
eigenmode solution, viscous eddies will appear for angles smaller than 79.6° (Moffatt 



1964); in that case the real part, TZ(Xe) > 3.8, so that the eigenmode is subdominant 
with respect to the inhomogeneous solutions. Hence, in the evaporation problem, the 
eigenmodes become important only for angles where no eddies are present. 



3.3. Lubrication limit: 9 <C 1 

Now we have the full analytical solution to the evaporation-driven velocity field in a 
wedge available, we can check whether for sufficiently flat droplets (i.e. small contact 
angles) the lubrication approximation can be applied. Therefore, we expand the Stokes 
flow solution for the flux condition for small 9, <f>, yielding 



u p {p,<l>)~3A(6)U 



1 



1 



9 (p/Ry-W) 



(3.3) 



In the lubrication approximation, we can write the radial velocity in terms of the radial 
distance from the drop center, r, distance to the solid substrate, z, and drop height h, 
as flHu fc Larsorl||2005| |Marin et aZ."||2011| > 

3\/2, 



u r (r,z) 



-U- 



v/1 - (r/R) 



(3.4) 



For small enough 9, A{9) — > V^/n |Popov 2005) and p ~ 1 — r/R. Furthermore, since 
z = psin0 ~ pcf) and h = psin9 ~ pO, z/h ~ <j>/0, and (3.3) leads to the same expression 



as (3.4). Hence, the Stokes solution converges to the lubrication solution for small 9. To 



illustrate this, we plotted the lubrication solution together with the Stokes solution for 
different contact angles; see figure [7j For convenience, we normalized all velocities by the 
velocity at the interface. From these results it is clear that even close to the contact line, 
where the evaporative flux is diverging, the lubrication approximation can be applied. 
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Figure 7. (a) (Colour online) The radial velocity u g ( </>), normalized by the velocity at the 
free surfac e u p (6), in the lubrication approximation (3.4 1 (solid black line), and the Stokes flow 
solutions ( ED) (dashed lines) for 6 = 20° (blue), £T=30° (red), 9 = 40° (green), 6 = 50° 
(orange), (b) Streamline plot of the lubrication solution for 9 = 50°. The velocity in azimuthal 
direction is derived from the continuity equation. 



4. Discussion 

We have derived the analytical solution to the Stokes-flow problem in a wedge geom- 
etry with an evaporative flux boundary condition. From the Stokes-flow solution it was 
found that the lubrication approximation accurately describes the velocity field in an 
evaporating drop all the way down to the contact line, for small enough contact angle. 
For larger contact angles, reversing flow structures are observed, and the eigenmode so- 
lution of the problem becomes of importance. Hence, the direction of the flow (towards 
or away from the contact line) depends on the contact angle of the droplet. The contact 
angles at which new separatrices in the velocity field appear are calculated from the 
exact solutions. This flow reversal, the appearance of separatrices and the dominance 
of the eigenmode solution for evaporation-driven flow have not been reported before in 
numerical studies ( Fischer||2002 Hu & Larson 2005), which focussed on smaller angles, 
nor in analytical studies (Masoud & Felskc 2009 Petsi & Burganos 2008), where differ- 
ent boundary conditions were used. It would be interesting to further investigate these 
aspects in numerical simulations of the flow in the entire drop, in particular near the 
critical point, where the exponents of all three contributions (flux, hinge and eigenmode) 
to the total solution are very close. 

The direction of the surface flow in the evaporating drops changes when the contact 
angle passes through the critical angle 9 C . In addition, Marangoni stresses induced by the 
non-uniform evaporative flux can change the surface flow direction. We can re-examine 
the influence of the Marangoni effect by direct comparison of our results to the solution 



with a Marangoni boundary condition, presented by Ristenpart et al. (20071. In our case, 



hence in absence of Marangoni stresses, the radial velocity scales as 



whereas the solution for the Marangoni case scales as ( Ristenpart et aL]|2007 ) 



? 2 [/Ma/5 A , 



(4.1) 



(4.2) 



where Ma= pi(3RAH v / \ik the Marangoni number, with AH V the specific latent heat of 
evaporation, /? the dependence of surface tension on the temperature, and k the thermal 
conductivity of the liquid. For water, the Marangoni number is of order 10 5 . The factors 
1/6 and 2 only appear in the small contact-angle limit, otherwise, the ^-dependent 



prefactors are of order unity. From (4.1 ) and (4.2 ) it is clear that the two expressions have 



a different power in p, and therefore a cross-over length scale exists. Below this critical 
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length p Cl the lowest power, hence the evaporation-driven solution (4.1 1, will dominate, 
above this length the Marangoni-driven flow is dominant. The cross-over length scale is 
given by 

hence, this length scale is of order 10 _8 m for a millimetre-sized droplet, as long as the 
contact angle is of order 1. Only for extremely small contact angles, the length scale over 
which evaporation dominates the Marangoni effect becomes comparable to the drop size. 
Hence, the Marangoni-effect should alter the velocity pattern in evaporating droplets 
significantly, as described by the solutions obtained by Ristenpart et al. (2007). The 
question why this is not observed experimentally for water drops remains, with as possible 
explanation that it has to do with surface- active contamination which could suppress the 
Marangoni effect ( |Hu k Larson][2006| . 

We described flow structures inside drops on partially wetting substrates, in cases 
where the contact line is pinned. When the contact line is not pinned but free to move, 
this will lead to an additional velocity field inside the drop, as described by |Huh fc| 



Scriven (1971). Since the Stokes flow problem is linear, the moving contact-line solution 



by Huh k Scriven ( 1971 ) could just be superimposed on our evaporation-driven solution. 
In that case, two velocity fields that are in opposite direction will arise in the problem: 
one with a mean flow directed towards the contact line (for contact angles below 127°), 
driven by the evaporative flux from the drop surface, and one with a mean flow away 
from the contact line driven by the receding contact line. Since the evaporation-driven 
velocity field diverges close to the contact line for 8 < 90°, it will dominate the solution 
close to the contact line. A cross-over length can be defined where both velocity fields 
cancel each other (Berteloot et al. 2008| ). This length depends on the contact angle and 
can be as large as 100 ^tm: for small contact angles the entire flow field in the wedge close 
to the contact line will be dominated by evaporation. 

Interestingly, the diverging evaporative flux not only leads to a diverging velocity 
field inside the droplet, but also to a diverging pressure field, which scales as p ~ 
l/e(r]U/R)p x - 2 . For water drops, with U = 10" 6 m/s, rj = 10~ 3 Pa s, R = 10~ 3 m, this 
gives on very small length scales p = 10~ 6 (p = 10~ 9 m) that p = l0 3 /8 Pa. This pressure 
is large compared to the Laplace pressure p 1 = -f/R = 10 2 Pa. Therefore, we expect on 
small scales (p ~ nm) some additional curvature of the interface. This subtle problem 
has been studied in detail by Berteloot et al. (2008): for angles ~ 0.1 a change in contact 
angle of about 15% was reported. For angles of order unity, however, no significant effect 
of the evaporation-driven flow on the contact angle was found. Note that, since p ~ p x ~ 2 , 
the pressure will diverge for all contact angles, since the leading order exponent A < 2 
(see figure [2J). For 9 < 90°, or A < 1 , the evaporation prob lem is even more singular than 
the moving contact-line problem by Huh k Scriven ( 1971 ), for which p ~ p^ 1 . Moreover, 
the introduction of a slip length does not change the pressure exponent given by the 
evaporativc-flux boundary condition, contrarily to the case of a moving contact line of a 
non-volatile liquid. 

Finally, it is still an open question what regularizes the radially diverging velocity in- 
side an evaporating drop with a contact angle smaller than 90°. In previous studies, two 
approaches have been used. Either the evaporative flux singularity is smoothed mathe- 



matically (Fischer 2002 Masoud k Felske 2009 Petsi k Burganos 2008 Poulard et al. 



2005), or by the introduction of a precursor film, which regularizes the singular wedge 
geometry near the contact line ( Eggers k Pismen|2010 Pham et a/.||2010 Poulard et al. 



2005 Semenov et al. 2011 1. However, for droplets on a hydrophobic substrate, a precursor 
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film is not expected, and a proper explanation for the regularization of the evaporative 
flux is still lacking. The large pressure close to the contact line could reduce the vapour 
concentration, according to the Kelvin equation ( Eggers fc Pismen|2010 ) but for contact 



angles of order unity the pressure appears insufficient for this effect, even at the nanome- 
tre scale. Furthermore, if the rate at which molecules are transported by diffusion exceeds 
the rate at which they transfer across the interface, the air just above the droplet is no 
longer saturated with vapour and the evaporative flux no longer diverges. We speculate 
that another option could be that the regularization occurs on the length scale of the 
mean free path of the water vapour molecules in the surrounding air, which is of the 
order of 100 nm, and sets a lower bound on the validity of the diffusion problem. 

We are grateful to Alvaro G. Marin, Detlef Lohse, Pierre Colinet and Howard Stone 
for valuable discussions. We acknowledge the financial support of the NWO-Spinoza 
program. 
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